N90- 10386 


THE 
I DENT IF 


GENERALIZED DRIFT FLUX APPROACH: 
CATION OF THE VOID-DRIFT CLOSURE LAW 


j. A. Boure 


Commissariat a l'Energie Atomique 
Centre d'Etudes Nucleaires de Grenoble 
Service d'Etudes Thermohydraul iques 
85 X, F 380^1, Grenoble Cedex, France 


abstract 


* “ al " »nd the potential advantages of 

generalized drift flux models are recalled. In particular it is 
stressed that the issue on the propagation properties and on the 
mathematical nature (hyperbolic or not) of the model and the 
problem of closure are easier to tackle than in two- fluid models. 

The problem of identifying the differential void-drift closure law 

inherent to generalized drift flux models is then addressed. Such 

a void-drift closure, based on wave properties, is proposed for 

u bly flows. It involves a drift relaxation time which is of che 
order of 0.25 s. 


is observed that, although wave properties provide essential 
closure validity tests, they do not represent an easily usable 
source of quantitative information on the closure laws. 
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I. introduction 


Generalized drift flux models were recently shown (Boure, 1988b) 
to be attractive alternatives for current 1-D two. fluid models . 
Drift flux models are characterized by the uses of a single 
momentum balance (the mixture balance instead of two phasic 
balances) and of a void-drift closure law. In classical drift flux 
models the void , drift closure law is expressed through an 
algebraic equation, which amounts to ignoring nonequilibrium drift 
effects. In generalized drift flux models, the void-drift closure 
equation is a partial differential equation. 

Generalized drift flux models and two-fluid models are compared in 
the next section and in table 1. The comparison brings out the 
drawbacks of two-fluid models. However, both kinds of model must 
be complemented by closure laws. In particular, generalized drift 
flux models need a void . drift closure law which remains to be 

specified . 

Since generalized drift flux models were introduced to account for 
the properties of kinematic waves (Boure, 1988a), it seems 
logical to use these properties to identify (i.e. to evaluate the 
coefficients of) the void-drift closure equation. The purpose of 
the present paper is to discuss the identification problem, 
assuming that the void-drift closure equation may be approached by 
a quasi linear differential equation of the first order. 

H. a rf.MTNDER on GENERA L IZED DRIFT FLUX MODELS AND, TWO-FLUI D 
MODELS 


The comparison between generalized drift flux models and two 
models is summarized in table 1 which, like must of the substance 
of this section, is taken from Boure (1988b). The mass and energy 
balances are parts of both kinds of models, and they are not 
discussed further hereafter. It is only noted for completeness 
that they require a few closure equations, in particular for the 
mass and energy transfers at the walls and at the interfaces. 
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The momentum balances contain the two phasic averaged pressures P Q 
and P L , the subscripts G and L corresponding to the two phases. 
For the following discussion, it is convenient to express P Q and 

P L in terms of the average pressure P and the pressure difference 
P LG , with : 


P = a P G + (1 - a) P L , 


P LG = P G - P L 


( 2 . 1 ) 


a being the void fraction. Now, in the set of balance equations, 
the two phasic momentum balances are equivalent to a subset of two 
equations, namely the mixture momentum balance and the pressure 

difference equation, obtained on eliminating P between the two 
phasic balances. 


The mixture momentum balance may be written : 


dP (* w G 

dz + “ P G [JT 


d W, 


+ W„ 


G v "G ) (d w a W. \ 

at + — J + d - “> Pt (— + w t — ) 


aR 3r 

+ aT + dz + MlG (W G " W L ) _ I cr + F w + pg = 0 


( 2 . 2 ) 


the terms of the second line representing respectively : 

- a term accounting for fluctuation and transverse distribution 
effects 

- a term accounting for longitudinal stress variations 

- an interfacial mass transfer term 

- a surface tension term (interfacial) 

- a friction term (walls) 

- a gravity term. 

t and z are respectively the time and the space variables, p G and 
P L are the phase densities, W, and W, are the phase average 
velocities, g is the gravity acceleration, and : 

p = a p G + ( 1 - a) p L (2.3) 
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w being the local instantaneous velocity along Oz and j the local 

k ... 

instantaneous deviatoric stress tensor, indicating the 

conditional time or ensemble averaging operator and < > the 

space averaging operator, and n being the unit vector of the Oz 
axis, R and t are defined as : 


R 


-G L 

= (a. p G (w Q - W G ) 2 + (1 - oc) p L (w L - W L ) 2 ^ (2.4) 


T = 


/ -G -L 

- ^ a x + (1 “ a ) I, 


"L 


n ) • n 
z / z 


(2.5) 


Independently of the mass transfer, already present in the mass 
balances, eq. 2.2 requires four closure laws for R, t, 1^ and Fy . 
Turbulence effects are present through R. 


The pressure difference equation may be written : 
W, 


a (1 -a) 


p G [' 
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+ w r — ~ 
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- Pl 

( d W L 

■ + w L 

3 W L x 3 P LG 


i at 
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' W GI> 

+ a 

(W L - 

W LI>] 

(2.6) 


+ L, + L fu + L F| - a (1 - a) p GL g = 0 


the terms of the second and third lines representing respectively 

- a term accounting for fluctuation and transverse distribution 
effects 

- a term accounting for longitudinal stress variations 

- an interfacial mass transfer term (W Q , and W L , are interfacial 
averages of the phasic velocities) 

- an induced inertia ("added mass") term 

- two friction terms (respectively wall and interfaces) 

- a gravity term, with 


Pgl = Pl - p G 


(2.7) 


The above terms are defined in Boure (1988b). For instance : 


o 


(2.8) 


** - (1_a) 57 (“ Pg < W G _ U G> 2 ) “ “ 57 (< 1_ 0l> P|. ("l" \ > 2 ) 

L, 4 - (l-o.) f(oc /■ n_) • n ; + a ± ((1-a) / • „/ . „ (2.9) 

Equation (2.6) requires seven closure laws for T t i.t 

T? ' T 1 > "g i > 

W ll- L i - l fw * L f i ■ Moreover, retaining eq. (2.6) implies, at least 
in theory, retaining also the interfacial momentum balance, which 
requires one supplementary closure law. 

In two-fluid models, eqs (2.2) and (2.6) are both implicitly 
written. A crucial point is that eq. (2.6) is stiff. This means 
that, due to the relative orders of magnitude of its terms, small 
variations on the closure laws (and in particular on the induced 
inertia and interfacial friction closure laws) induce large 
variations on P LQ (when P LQ is not merely ignored) and/or W Q . 
Hence the difficulty to adjust the very badly known closure laws 
of eq. (2.6) to avoid unrealistic values of the light phase 
velocity . 

A second crucial point is that, whenever eq. (2.6) is closed by 
algebraic laws only, a non-hyperbolic set results : two 
characteristic velocities are complex conjugate, with the 
consequence that the model is unconditionally unstable. 

In drift flux models, eq . (2.6) is not written, which is 
acceptable since the value of P LG does not really matter in 
practice. The two foregoing difficulties are not encountered. 

Besides the closure laws already mentioned, current two -fluid 
models are closed through an assumption on P^ (the set of 
6 phasic balance equations involves 7 dependent variables, namely 
two pressures, two velocities, two enthalpies and the void 
fraction) . Such an assumption imposes an artificial constraint on 
the pressures and pressure gradients. It disturbs the description 
of the corresponding propagation phenomena. 
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Generalized drift flux models are closed through the direct 
description of the void- drift dynamic dependency. Assuming a 
quasi-linear , partial differential relationship of the first 
order, and using the convenient variables : 

W * a W G + (1 - a) W L , S £ a (1 - a) (W„ - W L ) = a (W c - W) (2.10) 


(center of volume velocity and drift), it may be approached by 
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( 2 . 11 ) 


Equation (2.11) requires seven closure laws, i.e. the same number 
as eq. (2.6) but with simpler physical significances and more 
straightforward consequences : f is the fully-developed drift 
value, 9 a relaxation time, I and U are respectively the sum and 
product of the two characteristic velocities corresponding to the 
kinematic waves, t] expresses inertia effects, W 2 and W 4 are 
averaged velocities close to W. 

It can be concluded that developing a closure set for eq . (2.11) 
and using generalized drift flux models appear as less hazardous 
and hopefully easier than developing a closure set for eq. (2.6) 
and using current two- fluid models. 


m TWVT.ITF.NCE OF THE VOID - HR TFT CLOSURE ON THE PROPERTIES OF 
VTNF.MATTC WAVES ("Direct problem") 


As long as the kinematic wave velocities are small with respect to 
the sonic velocity, the properties of small harmonic kinematic 

waves of the form 
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X 


i (cot - kz) 


(3.1) 


e 


where x Q (constant) (0 and k are real or complex quantities, result 
from the approximate dispersion equation (Boure, 1988b) 

w " C a k + i 0 (go 2 - Ztok + IIk 2 ) = 0 (3.2) 


where : 


C 


a 


W + 



(3.3) 


A first consequence is that the kinematic wave properties, which 
do not significantly depend on rj , W 2 , , cannot be used to 
evaluate these quantities (r| is related to the sonic velocity, W 2 
and W 4 have only weak influences) . 

In the exploitable experiments (Tournaire, 1987, Boure, 1988a) co 
is imposed and the kinematic wave velocities V and their spatial 
amplification coefficients k ? result from the data processing. 
Eqs . (3.1) and (3.2) must therefore be used with co real and : 

k = k r + i kj (3.4) 


from which 


k { z i (cot - k z) 

x=x 0 e'e (3.5) 


co 



k p , which does not depend on the frame of reference, can be used 
instead of co, which does, to characterize a wave. 


Introducing eqs. (3.4) and (3.6) in the dispersion equation (3.2) 
and separating the real and imaginary parts lead to : 
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k r [V - C, + 9 S V kj - 2 8 n kj ] = 0 


(3.7) 


- C a kj + 9 V 2 kj; - 0 I V kj: + 0 n (kj: - k?) = 0 (3.8) 

Equations (3.7) and (3.8) provide the solutions to the direct 
problem, viz. computing the properties of the kinematic waves of 
wave number k p when C a , 9, Z and II are known. There are two 
solutions corresponding to two modes (noted with the subscripts 3 
and 4) . In the experiments it was found that modes 3 and 4 are 
respectively predominant at low void fractions (0 < a < 0.25) and 
at "large" void fractions (a > 0.30) 


In particular for k p = 0 , the two solutions are : 


k- =0 with — 


kj 9 <C a - C 3 ) (C* - C 4 ) 


V = C„ 


(3.9) 


and 


k i = 


e n ' v Ca 


i^i j. i_ 

C ; c 3 C 4 c a 


(3.10) 


Cj and C 4 being defined by : 


2 & C 3 + C 4 


n a c 3 c 4 


(3.11) 


For k -t> oo, the two solutions are : 
r ’ 


V - c 3 k, = - 


- C 3 


0 Cj (C 4 Cj ) 


V - c 4 kj = - 


- c a 


004(^4 ^3 ) 


(3.12) 


(3.13) 
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In equations (3.7) and (3.8) the three quantities 0 k p , are 
present only through the two products (0 k p ) and (0 k. ) . In actual 
experimental runs, k p is never zero and the equations may be 
written : 

V - C 0C + (C 3 + c 4 ) v (0 k i ) - 2 C 3 C 4 (0 k f ) 3 0 (3.14) 

9 k i (Q k, ) 2 

- c a + V 2 - (Cj + c 4 ) V + c 3 c 4 - c 3 c 4 = 0 (3.15) 

(6 k p ) (0 k r ) 2 


Since in the exploitable experimental data, V does not 
significantly depend on k r (no' significant dispersion), it is 
convenient to eliminate V between eqs . (3.14) and (3.15). Eq . 
(3.14) yields : 

C a + 2 C 3 C 4 (0 k- ) 

V = (3.16) 

1 + (C 3 + c 4 ) (0 k f ) 



The solutions for 0 k- are then the real solutions of the 
equation 


2 C a - (C 3 + C 4 ) - <C 4 - C 3 ) 2 (0 k, ) 
1 + (C 3 + C 4 ) 0 k, 



Computing directly 0k* as a function of 9 k p from eq. (3.18) is 
not straightforward. It is more convenient to transform eq. (3.18) 
to express 0 k p as a function of 0 k. 
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(3.19) 


(0 k r ) 2 
0 k , 


[c a + c 3 c 4 e k, ] [i + (c 3 + c 4 > e k, ] 2 

( C a- C 3> ( C a — C 4> " ( C 4~ k i l C <* + C 3 C 4 0 k i 1 


For a given value of k- , eq. (3.19) yields zero or one real 
positive value of k p . 


IV. IDENTIFICATION OF THE VOID-DRIFT CLOSURE EQUATION FROM 
KINEMATIC WAVE PROPERTIES ("Inverse problem") 


The problem posed in this paper is the determination of C a , 9, Z 
and n, using the experimental data on kinematic waves. It is the 
inverse of the problem of section 3. 


The principle of the method is to write eqs. (3.7) and (3.8), for 

instance, for several sets of experimental conditions for which 

k V and k- are known and to use the resulting equations to 
r * t 

compute C a , 9, Z and IT. 

In the exploitable data, as already noted, V does not depend 
significantly on k p . Accordingly, when a single mode is 

predominant, V may be expected to be close to both C 3 (or ) and 
C a (or C^) . This has two consequences : 

1. Whenever a single mode is predominant, the experimental 

correlation for V should be a correlation for C a as well. This is 
corroborated by the fact that is satisfies the definition (3.3). 
In the experimental conditions of Tournaire (1987) (upward 
vertical flow, low pressure), it leads to (Boure, 1988a, f and 
C a - W in m/s) : 


For 0 < a < 0.2 (mode 3 predominant) 

f = 0.22 a (1 - a) [1 - 1.25 a (1 - a)] 

C a - W = 0.22 (1 - 2a) [1 - 2.5 a (1 - a)] 


(4.1) 
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For 0.3 < a < 0.41 (mode 4 predominant) 


f = 0.22 a - 0.028 

c a - W = 0.22 


(4.2) 


For 0.2 < a < 0.3 (the two modes coexist) the values of f and 
C a - W may be interpolated between (4.1) and (4.2) with, from the 
experimental data : 


for a = 0.25 


C a - W s- 0.08 


(4.3) 


2. When mode 3 (respectively mode 4) is predominant, C a - C 3 
(respectively C 4 - C a ) should be "small". 

being known and V eliminated, the problem may now be 
reformulated, eqs . (3.7) and (3.8) being replaced by eq. (3.18) or 
(3.19) to be solved for 8, C 3 , C 4 . Only mode 3 results are 
exploitable since mode 4 results for kj are too few in number. In 
view of the forms of eqs. (3.18) and (3.19), the foregoing problem 
is far from simple. This is confirmed by fig. 1 in which, as 
suggested by eq. (3.19), the experimental results for - kj|/k. are 
plotted as a function of — k- , and which exhibits an important 
scatter (in the representation of fig. 1, the points corresponding 
to Ik,. I < 0.1 are subject to large errors and therefore 
meaningless. They were not plotted in the figure). 

In the domain in which mode 3 is predominant, the conditions : 

C 3 < C a < C ot < C 4 (4.4) 

may be expected to hold (see Boure , 1988a). They entail 
C °c - C 3 ^ 1 ^ - c a c a 

0 c 3 (C 4 - C 3 ) e (c 3 + c 4 ) 0 c 4 (C 4 - C 3 ) < 6 C 3 c 4 (4 ' 5) 
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Then, as k f decreases from zero, kjl/kj resulting from eq . (3.19) 
starts from the value corresponding to eq . (3.9) with a slope 
which may be of either sign. However, when Ikj I is sufficient, 
k^/k- also decreases. It tends towards - oo for the value 
corresponding to eq. (3.12). For the values of k- comprised 
between those given by eqs . (3.12) and (3.13), there is no 
physical solution (k 2 r < o) . Finally, for the values of k f 
comprised between those given by eqs. (3.13) and (3.10), there is 
a solution again corresponding to mode 4. 

By trial and error, 0, and C 4 may be adjusted to fit the curve 
representing eq . (3.19) to mode 3 data. In view of the absence of 
dispersion C 5 may be expected to be close to C a . 

For 0 < a < 0.20 : the following set of values is acceptable 

0 - 0.25 s ) 

C ot C 3 = °* 02 V s < 4 * 6) 

C 4 - C a = 0.08 m/sJ 

Precise adjustement would need more accurate data for the wave 
velocities and especially for the damp ing/amplif icat ion 
coefficients. Such data does not exist and cannot be expected to 
be obtained soon in view of the available instrumentation. For 
a > 0.25, no sufficiently accurate data is available. 

Equations (4.1) and (4.6) confirm that, even at low pressure, the 
relevant velocity differences corresponding to fully developed 
conditions are fairly small in bubbly flows. They are negligible 
as soon as W is large enough (say 3 m/s). Accordingly the 
correlations for C a - W, C 3 - W, C 4 - W are probably not crucial. 
On the other hand, the drift relaxation time 0 is an essential 
parameter of the generalized drift flux model. 
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V. CONCLUSIONS 


After a reminder on the main characteristics and the potential 
advantages of generalized drift flux models, the problem of the 
identification of the differential void-drift closure law they 
imply has been addressed. 

Two advantages of generalized drift flux models versus complete 
two- fluid models are : 

1. The correct description of the kinematic wave phenomena and the 
straightforward control and interpretation of the mathematical 
nature (hyperbolic or not) of the model set of partial 
differential equations that they enable. 

2. The simplification of the closure problem, involving only 
closure laws of simple physical significance and easy to assess. 

On the other hand, it has been found that the available data on 
kinematic wave properties is not quite adequate to enable the 
identification of the void-drift closure law. Such an 
identification would need a very good accuracy (difficult to reach 
in practice) on the wave damping or amplification coefficients. 
Accordingly kinematic wave properties seem to be more useful as a 
closure validity test than as a source of quantitative information 
on the closure laws. 

However, a void- drift closure, suitable for bubbly flows, has been 
adjusted on the available kinematic wave data. It involves several 
velocities which differ only slightly from each other and from the 
average fluid velocity but whichare necessary to the description 
of the kinematic wave properties. It also involves a drift 
relaxation time which is an essential parameter of generalized 
drift flux models and which was tentatively found to be of the 
order of 0.25 s for the upward flow of air-water mixtures at low 
pressure . 
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Table 1 

Simplified comparison of modeling strategies 


Drift flux approach 


Common features 


Two-fluid approach 
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